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based on the classic DTC working with a constant switching frequency of the inverter. 
Another remarkable aspect is the complexity of the method proposed, both in the 
control unit of the inverter and in the number of correctors necessary for the control 
of the torque. The selection table and hysteresis have been eliminated. This method 
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significantly improves the torque and flux oscillations and improves the dynamics of 
Direct-backstepping control the drive by making it less sensitive to load torque disturbances. The proposed method 
Direct torque control is developed and designed using Matlab/SIMULINK to show the effectiveness and 
Multiphase motor performances of the DTC-Backstepping. 
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1. INTRODUCTION 

Multiphase machines offer an interesting alternative to the reduction of the stresses applied to the 
switches as to the windings. Indeed, the multiplication of the number of phases allows a fractionation of 
the power and hence a reduction of the voltages switched to a given current. Moreover, these machines 
reduce the amplitude and increase the frequency of the torque ripples, thus enabling the mechanical load to 
filter them more easily. Finally, the multiplication of the number of phases offers increased reliability by 
allowing operation when one or more phases are open. Multiphase machines are present in the fields of marine, 
railway traction, petrochemical industry, avionics, automobile, etc. The rotor flux oriented control mostly 
employed in multiphase motor drives is based on effective control of the magnetic state. However, this structure 
generally requires the installation of a speed/position sensor on the shaft. Moreover, it remains very sensitive 
to the variations of the parameters of the machine, many research have developed this strategy of control 
applied for multiphase machines [1-6]. Direct torque control (DTC) has been initially developed for induction 
machines in the years 1986 and 1988 by TAKAHASHI and DEPENBROCK. In the literature there is an 
extensive research on DTC applied for multiphase machines [7-14]. In [7,12] a DTC strategy for dual 
three-phase induction motor drives is discussed. In [8] describes an investigation of direct torque control 
for multi-phase permanent magnet synchronous motor drives. A DTC for five-phase synchronous motor has 
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been developed in [9]. Also in [11], authors established a DTC algorithms for a split-phase induction motor. 
Compared to vector control, the DTC control is much less sensitive to parametric variations and allows to ob- 
tain dynamics faster control. A good accuracy in measuring the position of the rotor is not necessary since 
only the sector in which the flux is located matters to determine the configuration to be used.On the other 
hand, the flux and torque of the machine must be estimated or observed using different ways. The synthesis of 
such estimators is not trivial and constitutes a difficulty for the implementation of this control [15]. By nature, 
torque oscillations exist in DTC drive. The reduction of the hysteresis bands with a given sampling period does 
not always have an effect on the amplitude of the torque oscillations. In this case, in order to reduce torque 
oscillations, it is necessary to decrease the sampling period. The switching frequency is not controlled, it varies 
according to the operating point. DTC with constant switching frequency using SVPWM is reported in [16,17], 
however, the algorithm is more complex, nevertheless, torque and flux oscillations are reduced. One find in [18], 
the extended DTC, this control uses selection table of optimal vectors and considers an additional input: 
the sign of the change of electromagnetic torque. In the extended DTC, the torque regulator has a three- 
level output, unlike the conventional DTC where only two levels were considered. One can observe a reduction 
in torque oscillations. In [19,20], the authors present an algorithm allowing to have a constant switching 
frequency, its main feature is the removal of the hysteresis regulators and the vector selection table, which 
eliminates the problems associated with them. With this control method, the inverter operates at a constant 
frequency, since a PWM technique is applied. The objective of this method is to realize a direct control of 
the stator flux vector in a reference frame linked to the stator, the polar components of these two vectors are 
obtained by their projections on the reference frame. From these components, the desired stator flux vector at a 
given instant is calculated. The PWM is applied to this vector to obtain the switching states of the inverter. 
Applications using fuzzy systems have been developed in [21,22], where the hysteresis blocks have been 
replaced by fuzzy controllers.vFuzzy logic is used to achieve a compromise between torque control and flux 
control, and they do not require an exact mathematical model of the machine. However, it draws strong current 
pulses which normally translates into higher torque ripple. The sliding mode control (SMC), due to its robust- 
ness to uncertainties and external perturbations, can be applied to uncertain and disturbed nonlinear systems. 
It is a question of defining a so-called sliding surface according to the states of the system. The synthesized 
global command consists of two terms: the first allows to approach this surface, the second allows the reten- 
tion and the sliding along it [22]. The combination of SMC and DTC control reduces torque and flux ripples. 
Its main feature is the removal of the hysteresis regulators and the switching table, which eliminates the prob- 
lems associated with them. The disadvantage of this association (DTC-SMG)is the use of the saturation func- 
tion that introduces a static error that persists as well as the need to have a knowledge of the dynamics of the 
system. In [23], authors addresse a new approach to adapt the concept of the Predictive Control Model(MPC) to 
Direct Torque Control (DTC) in the five-phase induction motor control. The proposed algorithm improves the 
performance of a DTC controller by retaining the electromagnetic torque and the stator flux modulus within 
predefined hysteresis bands while minimizing switching losses. The MPC controller extracts the strings of 
switching sequences on the forecast horizon. Dynamic programming is implemented to select the switching 
sequences that minimize the cost function on power losses. Both simulations and experimental are carried 
out on specific models and the results verify the advantages of the proposed DTC method in comparison with 
conventional DTC. Model-based predictive direct control methods are advanced control strategies in the field 
of power electronics to control induction machines [24]. The torque control prediction (PTC) method evaluates 
the electromagnetic torque and the stator flux in the cost function. The switching vector selected for use in insu- 
lated gate bipolar transistors (IGBT) which minimizes error between references and predicted values. System 
constraints can be easily understood.Behavior and robustness and transient performance are evaluated. Due to 
the inconvenient of the DTC control such as high torque and currents ripples push us to think and to design 
new control to improve the direct torque control technique. The idea in this paper is to use the combination 
between the DTC control and a non-linear control in this paper is the Backstepping control. As the authors 
know there is no study in the literature treated this approach applied to the multiphase motor and specially 
five-phase induction motor. 


2. DTC-BACKSTEPPING CONTROL OF FIVE-PHASE INDUCTION MOTOR 

The induction motor drive system used to develop the proposed controller strategy is composed mainly 
of five-phase voltage source inverter and a symmetrical five-phase induction motor with phase shift of the 
windings equal to 27/5. All the components of the drive system are schematically illustrated in the Figure 1. 


Torque ripples improvement of direct torque... (Hamdi Echeikh) 


66 o ISSN: 2088-8694 


Using various assumptions, such as negligible core losses, uniform air gap, sinusoidal MMF distribution and 
symmetrical distributed windings, then the five-phase induction motor model using the vector space decom- 
position method presented in [25]. The model represented in two orthogonal planes, where the first plane is 
involved in the fundamental torque production (a — 8 representing the fundamental component and contains the 
harmonics of order 10n + 1,n = 0,1, 2, 3...). The second plane is connected to the losses in the motor (x — y 
represented the supply harmonics of the order 10n + 3,n = 0,1, 2, 3...). Finally, the zero-sequence component 
is not considered due to the neutral point isolation isolation(z -axis representing the supply harmonics of order 
5n,n = 0,1, 2,3...). Using the selection of the stator currents in two planes a — 8 and x — y and the rotor 
fluxes in the plane a — £ as state variables, £1 = isa, 22 = 158, 73 = Yra, T4 = Wr GB, V5 = tse and TE = isy. 
The motor drive equations can be expressed in the following form 








£1 = C1£3 + C2WT4 — €3X1 + C41 





£2 = Cı T4 C2WT3 C3 T2 F C4U2 





T3 = —C5 £3 — WT4 + C6£1 

T4 = —C5T4 + wWT3 + C6T2 d) 
T5 = —C7L5 + C8V3 

Te = —C7Xe + C8V4 








w = Mı (£3£2 — £4£1) — M2T, — m3w 


With the coefficients given by 
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The inputs signals of system above are the applied stator voltages vı = Usa, U2 = Usß, U3 = Usa 
and vı = Vsy. The equations represent the electrical speed w and all the motor parameters, M is the mutual 
inductance, R, is the stator resistance, R, the rotor resistance, L, is the rotor inductance, Ls is the stator 
inductance and Ls is the stator leakage inductance. One presented a novel state variables to describe the 
novel model to developed the DTC-Backstepping controller, however, the separation of the model of five- 
phase induction motor is released into two different parts, mechanical part and electrical part, the new state 
variables are given as follows 


Ty = VratsB = VrBlsc = %3%Q — 4T 
by = 3(25 + x1) (2) 
Xy = T321 + T412 
Where T, is a virtual torque, the relation between the electromagnetic torque and the virtual torque is 
Tem = MıTy. Wy is virtual flux, the relation between the rotor flux and the virtual flux is Yw, = 2, where 


p2 = pa + Weng. In practice the components Yra and y,g are derived from the five-phase induction motor 
voltage model or current model given by 


sa = f sas i Rsisa)dt 


(3) 
Ysg = J (Usa — Rsisp)dt 
The rotor flux can be derived as following 
Yra = L, Ysa — OL ssa 
( ) (4) 


Ure = Lr (sp A oL sis) 


The five-phase induction motor model can be separated into electrical and mechanical parts using the new state 
variables given above, the mechanical part is 


w= mT, — mT, — m3w (5) 
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The derivative of the virtual torque is given by 
Ty = 1322 + 2329 — F421 — U4} (6) 
Replacing the derivative x1, x2, x3 and x4 from equation (1) one obtain: 
T,= — (c3 + c5 )(£3£2 — £41) — w(z3£1 + F442) — cow (x3 + £4) + c4(x3v2 — £401) (7) 
Using equation (2), the equation (7) can be rewritten as 
Ti, = —(¢3 + €5)Ty — wXy — 2c2w + caur (8) 
Where up = £302 — £401 = YraVsg — WreVsa The derivative of the virtual rotor flux Yẹ, gives: 
by = t3x3 + T44 (9) 
Using equation (1), equation (9) can be rewritten as 
Wy = —¢5(a2 + £2) + ce(£3£1 + £412) (10) 
Using equation (2), equation (10) becomes 
hy = —C5Xy + coo (11) 
The derivative of X, gives as follows 
Xy = F321 + £13 + Lyre + T42 (12) 
Using the equation (1), equation (12) becomes: 
X, = — (c3 + c5)(£3x£1 + 0402) + w(£3£2 — £421) + ce(? + £2) + ey (02 + £2) + c4(£3v1 + £402) (13) 
Replacing the expressions of equation (2) in equation (13), one obtains: 
Xo = —(c3 + ¢5)Xy + wT y + c6(2? + £3) + 2c1Yu + catey (14) 


Where uy = £301 +£4V2 = YraVı + Yrgv2 Indeed, for the mechanical part, the virtual torque T, is controlled 
with a virtual torque voltage ur and for the electrical part, the virtual rotor flux Yẹ, is controlled by the virtual 
rotor flux voltage uy. This will simplify the design of the DTC-backstepping controller in the following section 
based on the mechanical and electrical parts. One design by w* the electrical motor speed reference, the speed 
error is given by: 








zy=w*—w (15) 
The derivative of speed error can be rewritten as 
zy = w* -w= w* my Ty } mT, t m3W (16) 
Choosing the derivative of the error 2} = —k,z, where kı is a positive constant. One can get the reference of 
the virtual torque T, as 
T = (w* + mT, + mgw + kyz1)/m4 (17) 


The error of the virtual torque gives by z2 = T% — Ty, the derivative of z2 gives in equation (18). 





Zo = T* — T, = T* + (cg + e5)Ty twXy + 2equyy — caur (18) 
Choosing the derivative of the error z2 = —k22z2, one can get the virtual torque voltage ur as 
up = (T> + (c3 + ¢5)Ty + wXy + 2c2wpy + k2z2)/c4 (19) 


Using equation (17) the equation (19) can be rewritten as 


ur = ((w* + mT, + m3w — k?z1)/mı + (cg +¢5)Ty + wXy + 2cowy + k2z2)/c4 (20) 


Torque ripples improvement of direct torque... (Hamdi Echeikh) 


68 o ISSN: 2088-8694 


Finally, replacing the expression of w from equation (1), the expression of the virtual voltage ur 
given as follows 


ur = ((w* Line hy = m3MoT, — maw — kz +my1(m3+c3 +¢5)Ty)/my+wXy + 2cqwy +k22z2)/e4 (21) 


One assume that the virtual rotor flux of the five-phase induction motor given by Y , then its error is given by: 


z3 = Py — Yo (22) 
The dynamic of the error of the virtual rotor flux is given by the following equation: 
23 = i — Wo = Uh + csv — coXy (23) 
Choosing the dynamic of the error 23 = —k3z3, one obtains the virtual control X% as: 
Xi = (Pb + 2esthy + kazs)/co (24) 


Let is now define the error of the virtual control X, as: 


z4 =X- Xy (25) 
Its dynamic error is given by: l , 
Zi = X*— Xy (26) 
Using equation (14), equation (26) becomes 
ža = X* + (c3 + c5)Xv — WTy — ce(£? + £2) — 2cryy — cary (27) 
Choosing the dynamic of the error z4 = —k4za, where k4 is a positive constant, then the virtual rotor flux uy 
obtains as l 
uy = (XS + (c3 + c5)Xv — wT, = cele? + x2) — 2c1 py + k4z4)/c4 (28) 


Using equation (24), one can rewrite the equation (28) as 
Up = (be + 2c5%y — k3z3)/ce + (c3 + ¢5)Xy — wTy — cg(27 + £2) — 2ery + kaza)/c4 (29) 
Replacing the expression of Wy gives by equation (11) in equation (29), one obtains 
uy = (Hs +2c5(—2c5y 


From the two virtual voltages, both for torque and flux, one can acquire the stator voltages in the (a — 3) frame. 


c6 Xv) —k3z3)/c6 + (c3 +5) Xv —wWTy — ce(£? +25) — 2c1Yv + k4z4)/c4 (30) 





(31) 


U1 = Usa = (L3Uy = x3up)/2y 
V2 = Vsp = (sty + e3ur)/2py 


Now, One need to determinate the control voltages v3 and v4, then one define the following current errors 
25 and 26 


z5 = — lsr = 1 — T5 
i 9 Sx S£ (32) 


ae TAL teg 
26 = Usy — tsy = Usy — T6 

The dynamic of the currents errors zs and zę given as follow 
25 = thy ge = ee m e 
ee eae ae (33) 
26 = Usy — tsy = Usy — T6 


Choosing the derivative of the errors 75 = —k5z5 and zg = —kgze, the stator voltages in the reference (x — y) 
obtained in the following equation 


(34) 





fe = Usa = Ge. — C75 + k525)/cg 


V4 = Usy = (is, — C7 %6 T kez6)/cs 
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3. STABILITY ANALYSIS 

Usually it is difficult to find the asymptotic stability of time-varying systems because it is very 
difficult to find Lyapunov functions with a negative definite derivative. We know in case of autonomous 
(time invariant) systems, if V is negative semi-definite, then it is possible to know the asymptotic behaviors 
by invoking invariant-set theorems. However, the flexibility is not available for time-varying systems. This is 
where Babarlat’s lemma comes into picture. 


3.1. Lemma: Barbalt’s Lemma 
Suppose f(t) € C1(a,oo) and lim; f(t) = y where y < oo. If f’ is uniformly continuous, 
then lim: +. f’(t) = 0. 


3.2. Proof 

One will prove the result by contradiction. lim: f(t) # 0. Then > € > 0 and a monotone 
increasing sequence {tp} such that tn — co as n —> oo and | f’(t) |> € for all n € N. Since f'(t) 
is uniformly continuous, for such €, 5 x > 0 such that Vn € N. 


| t= tn |< x =| f'E) = f'n) IS €/2 


Hence if t € [tn, tn + x] then 
IFE = F n) = Fir) — FO) | 21 FO) | I F En) = FE) [2 e- €/2 = €/2 
Then since f(t) € C! one have 


tn +x ty tn+x tn tx tn +x 
f roa- | rodal f roa f rolaz f e/2dt = xe]? 


, a 
However 


tn +x tn 
im | f Oat- f° Od |= jim | Fn +) = Fl) | 


t-00 
=| lim f(t, +x) — Him fin) | =l71—-|7|= 0 
This is a contradiction. Therefore lim;,.. f’(t) = 0 


3.3. Application of the Lemma 
One assume the Lyapunov function is 


224 2b + 28 +23 4 22 + 22 








V = 35 
5 (35) 
The derivative of the Lyapunov function given by 
V = ZZ + 2229 + 2323 + Z4ž4 + 2525 + 2626 = kiz? kaz2 k3z3 kazi ksz5 kező (36) 
From equation (36), one have V < —kı eae then 
V(0) — V(co) > | —kızidt (37) 
0 
V is bounded. According to Barbalat’s lemma one can have 
lim z1 =0 (38) 
t— oo 
Also, one can have 
lim zo = lim z3 = lim z4 = lim z5 = lim z = 0 (39) 
t= oo t= oo t= oo t= oo t= oo 


At last, according equations (38) and (39), one can notice that motor speed , torque and rotor flux can track 
reference values asymptotically. 
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4. SIMULATIONS RESULTS AND DISCUSSION 

DTC-Backstepping control system block diagram is shown in Figure 1. Whenever the reference speed 
command w* is given, system compares it with the actual speed w. The error is used to determine the virtual 
reference torque T*. The reference virtual torque is compared with actual torque T. The torque error obtained 
is used to obtain the appropriate virtual voltage vector ur. The second reference input that is flux compare 
with actual flux. The error added to determinate the virtual control X. Like torque control loop, flux control 
loop is also operated to select the appropriate virtual voltage vector uy. However, After the determination of 
the two virtual voltage vectors which used to obtain the actual voltage vectors (a — (). The two other voltage 
vectors (x — y) are obtains from a small loop current controller. A program in Matlab/Simulink environment 
has been designed for both the voltage source inverter (VSI)-fed five-phase induction machine, and performed 
to analyze the viability of the developed control technique a number of simulations have been performed. 
A small prototype around 1.5 kW symmetrical five-phase induction machine with four poles has been used. 
Table 1 shows the parameters of the machine, which correspond to those of the real five-phase 
induction machine. 
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Figure 1. Scheme of Direct torque Backstepping Control of five-phase Induction Motor 


Table 1. Parameters of the five-phase Induction Machine 








Motor Data 
Parameters Symbol Value Units 
Stator resistance Rs 10 Q 
Rotor resistance Rr 6.3 Q 
Mutual Inductance M 0.42 H 
Stator leakage inductance Lis 0.04 H 
Rotor leakage inductance Lir 0.04 H 
Moment of inertia Jm 0.03 kg.m? 
Stator inductance Ls 0.46 H 
Rotor inductance Lp 0.46 H 
Number of pole pairs P 2 - 
Rated stator flux wn 1:27 Wb 
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Figures 2 and 3 summarize the simulation results. A start-up transient under no-load conditions, 
with speed references reversion from 400 rpm/min to -400 rpm, is considered, and the motor-speed, 
electromagnetic-torque, stator-flux, and stator-current responses and the torque ripple are shown for both 
methods (DTC and DTC-Backstepping). The obtained speed response using the DTC-Backstepping method 
is quicker than that using the DTC method, while the settling time is practically the same. By using DTC- 
Backstepping technique, torque ripples is reduced from 2.5 N.m to1.2 N.m as shown in figure.3. Hence torque 
ripple is reduced by 52%. Flux ripple is reduced from 0.1 Wb to 0.01 Wb as shown in figure.3. Therefore flux 
ripple is reduced by 90%. Itis observed that the DTC-Backstepping which based on non-linear approach allows 
effectively to reduce torque ripple as well as flux ripple. From the simulation figures, it can be concluded that 
torque and flux direct backstepping control has small torque, speed , flux and current ripple. Torque and flux 
direct backstepping control have better performance than classical direct torque control. The scheme can not 
only avoid the effect of rotor resistance to rotor flux orientation, but also reduce effect of hysteresis to torque 
and flux. From simulation results, one view that system design can achieve direct torque and flux control and 
motor speed tracking performance is better than classical direct torque control. The control scheme can realize 
system state variable decoupling efficiently. 
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Figure 2. Motor speed,torque, stator flux and currents left side DTC and right side DTC-backstepping 
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Figure 3. Stator flux revolution, torque and flux ripples, left side DTC and right side DTC-backstepping 


5. CONCLUSION 

This paper separates induction motor into mechanical part and electrical part. Backstepping coupling 
is applied to realize torque and flux decoupling. This can reduce the couple between state variable, and can 
make speed, torque and flux fast track reference values. System design is based on static coordinate of stator. 
The control scheme synthesize the direct torque and vector control, which not only realizes torque and rotor 
flux direct control, but also have the little torque and rotor flux ripple. Backstepping control design can not only 
provide better speed, torque and flux tracking performance, but also assert system robust performance under 
speed reversal. Moreover, the torque and flux are respectively reduced by 52% and 90%. 
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